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ABSTRACT 

We compare observations of AGB stars and predictions of the Elitzur & Ivezic (2001) steady- 
state radiatively driven dusty wind model. The model results are described by a set of similar- 
ity functions of a single independent variable, and imply general scaling relations among the 
system parameters. We find that the model properly reproduces various correlations among 
the observed quantities and demonstrate that dust drift through the gas has a major impact on 
the structure of most winds. From data for nearby oxygen-rich and carbon-rich mass-losing 
stars we find that (1) the dispersion in grain properties within each group is rather small; (2) 
both the dust cross-section per gas particle and the dust-to-gas mass ratio are similar for the 
two samples even though the stellar atmospheres and grain properties are very different; (3) 
the dust abundance in both outflows is significantly below the Galactic average, indicating 
that most of the Galactic dust is not Stardust — contrary to popular belief, but in support of 
Draine (2009). Our model results can be easily applied to recent massive data sets, such as the 
Spitzer SAGE survey of the Large Magellanic Cloud, and incorporated in galaxy evolution 
models. 

Key words: stars: AGB and post-AGB — stars: late-type — stars: winds, outflows — stars: 
mass-loss — dust, extinction — infrared: stars — 



1 INTRODUCTION 

Winds blown by stars on the asymptotic giant branch (AGB) are an 
important component of mass return into the interstellar medium 
and may account for a significant fraction of interstellar dust, in- 
cluding dust formation in the early universe (Sloan et al. 2009). 
Therefore in addition to its obvious significance for the theory 
of stellar evolution, the study of AGB winds has important im- 
plications for the structure and evolution of galaxies (Girardi & 
Marigo 2007; Marston et al. 2009). In our own Galaxy, its es- 
timated 200,000 AGB stars are a good tracer of dominant com- 
ponents, including the bulge (Whitelock & Feast 2000; lackson, 
Ivezic & Knapp 2002). AGB stars also have a great potential as 
distance indicators (Rejkuba 2004). 

Over the past few years, data of AGB stars have improved 
significantly. The 2MASS survey delivered an all-sky near-IR pho- 
tometric catalog to K < 14, and Spitzer obtained mid-IR (IRAC 
and MIPS) photometry for several thousand AGB stars from the 
Large Magellanic Cloud (LMC, Meixner et al. 2006). In addition, 
the MACHO survey provided high-quality light curves for ~ 22,000 
AGB stars from the LMC (Fraser et al. 2005), and the Northern 
Sky Variability Survey provided light curves for another ~9,000 
AGB stars brighter than V ~ 15.5 and with declination S > —38° 
(Wozniak et al. 2004). These new accurate and massive data sets 
are expected to rejuvenate studies of AGB stars, leading to model- 
based interpretation of diverse measurements including photome- 



try, outflow velocity, mass-loss rate, pulsations, etc. (e.g., Marigo 
et al. 2008). 

AGB stars are surrounded by an expanding envelope com- 
posed of gas and dust, that has a major impact on observed proper- 
ties. The complete description of such a dusty wind should start 
with a full dynamic atmosphere model and incorporate the pro- 
cesses that initiate the outflow and set the value of M. These pro- 
cesses are yet to be identified with certainty, the most promising 
are stellar pulsation (e.g. Bowen 1989; van Loon et al. 2008) and 
radiation pressure on the water molecules (e.g. Elitzur, Brown & 
Johnson 1989). Proper description of these processes should be fol- 
lowed by grain formation and growth, and subsequent wind dynam- 
ics. Two ambitious programs attempting to incorporate as many 
aspects of this formidable task as possible have been conducted 
over the last decade by groups at Berlin (see Winters et al. 2000; 
Wachter et al. 2008; and references therein) and Vienna (see Hofner 
1999; Dorfi et al. 2001; Hofner et al. 2003; Nowotny et al. 2005; 
and references therein). While much has been accomplished, the 
complexity of this undertaking necessitates simplifications such as 
a pulsating boundary. In spite of continuous progress, detailed un- 
derstanding of atmospheric dynamics and grain formation is still 
far from complete (e.g., Hofner & Andersen 2007). 

Fortunately, the full problem splits naturally to two parts, as 
recognized long ago by Goldreich & Scoville (1976). Once radi- 
ation pressure on the dust exceeds all other forces, the rapid ac- 
celeration to supersonic velocities decouples the outflow from the 
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earlier phases — the supersonic phase would be exactly the same 
in two different outflows if they had the same mass-loss rate and 
grain properties even if the grains were produced by entirely dif- 
ferent processes. Furthermore, these stages are controlled by pro- 
cesses that are much less dependent on detailed micro-physics, and 
are reasonably well understood. And since most observations probe 
only the supersonic phase, models devoted exclusively to this stage 
should reproduce the observable results while avoiding the pitfalls 
and uncertainties of dust formation and the wind initiation. 

In Elitzur & Ivezic (2001, hereafter paper I) we present the 
self-similarity solution of the dusty wind steady state supersonic 
phase. The model assumes steady-state spherically symmetric mass 
loss with prompt dust formation and no subsequent change in dust 
properties (i.e., no grain growth and sputtering). The included dy- 
namical effects are the radiation pressure force, dust drift through 
the gas, and the gravitational pull by the star. Paper I contains the 
description of the model and its solution in full mathematical rigor. 
Here we discuss observational implications of the model, and com- 
pare them to available data. The main questions that we ask are 

• Can the model explain the wind velocity profile, including the 
region of strong acceleration at small radii? 

• Does the implied density profile produce spectral energy dis- 
tributions in agreement with observations? 

• Do dynamical quantities, such as the wind terminal velocity 
and mass-loss rate, and photometric quantities such as colors and 
bolometric flux, display correlations predicted by the model? 

• How similar, or dissimilar, are dust properties inferred from 
observations of oxygen-rich and carbon-rich AGB stars? 

We first discuss dynamical quantities, and then focus on the 
spectral energy distribution. In §2 we describe the wind velocity 
profile. Section 3 presents correlations among final velocity Voo, 
mass loss rate M and luminosity L. In §4 we discuss the spectral 
energy distribution (SED) and its relationship with the wind dy- 
namical properties. In §5 we discuss the allowed range of physical 
parameters of dusty winds. The results are summarized and dis- 
cussed in §6. Equations from paper I are referred to by their number 
preceded by I. 



2 VELOCITY PROFILE 

The velocity profile for all winds, whether the dust is comprised 
of carbon or silicate grains, can be summarized with the simple 
analytic expression (see §3 in paper I) 
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Here 9 is angular distance from the star in the plane of the sky and 
do is a characteristic angular scale. The power k introduces a weak 
dependence on tv, the overall radial optical depth at visual. The 
small-Tv value k = | reflects the effect of the drift, which domi- 
nates the dynamics in that regime. At large rv, reddening effects a 
switch to a more moderate profile with k — 0.4. 

We discuss the wind terminal velocity in §3 below. Here 
we concentrate on the shape of the radial profile v/v x . The pro- 
file is mainly controlled by the characteristic angular scale 9o, cor- 
responding to the dust condensation radius. It can be determined 
from the measured bolometric flux, fboi, via 
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Table 1. Standard parameters for dust grains with size a = 0. 1 p.m: Qv is 
the efficiency factor for absorption at visual; Q* is the Planck average at 
the stellar temperature (2,500 K) of the efficiency coefficient for radiation 
pressure (see paper I, equation 4); \&n and m, defined in eq.fj] determine the 
quantity 4< that enters into the expression for the dust sublimation angular 
size (eq.fJJ- 



Here T c is the dust condensation temperature and $ is a di- 
mensionless function that sets the dust condensation radius as 
I r»* 1/,2 (r,/T c ) 2 , where r, and T* are, respectively, the stellar 
radius and effective temperature. All our model calculations as- 
sume T c = 700 K and a black-body star with T, = 2500 K. The 
function ^ is determined from the radiative transfer solution (see 
1.41 and 1.42) and can be adequately approximated by the simple 
analytic expression 



* = *o(1 + 0.005t v 1 ). 



(3) 



(2) 



The parameters *I/o and m are listed in Table Q] together with other 
relevant properties of the dust grains. In all the numerical calcula- 
tions here we employ amorphous carbon grains with optical proper- 
ties from Hanner (1988), and (the "warm" version of) silicate grains 
from Ossenkopf, Henning & Mathis (1992). We assume a single 
size a = 0.1 ^m to compute the grain absorption and scattering effi- 
ciencies; the grain size does not enter the solution of the dusty wind 
problem in any other way. From their effect on the variation of an- 
gular scale 6o with rv, the parameters $o and m introduce the only 
dependence of the velocity profile on dust properties. However, in 
practice this effect is rather weak: 9q varies by less than 10% for 
rv between and 20. 

We compare these model predictions with the spatially- 
resolved maser observations of three O-rich supergiants reported 
by Richards & Yates (1998). The combined measurements of SiO, 
H2O and OH masers trace the velocity profile over three orders 
of magnitude of radial distance from just outside the stellar at- 
mosphere. Comparison with the theoretical profile in equation Q] 
requires three input parameters. The terminal velocity is mea- 
sured directly in CO observations because most of the CO emission 
comes from the outer parts of the envelope, far beyond the acceler- 
ation zone. Next, we constrain the optical depth, tv, which controls 
the power-law index k, by fitting the spectral energy distributions, 
as described below in 3S](its effect on v(0) is minor). Finally, we 
determine the angular scale 8q from eq.|2]and the measured -Fboi. 
In using this relation we assumed T c = 700 K for all three stars, but 
the entire range 600-800 K is consistent with the data. All relevant 
parameters for the three stars are listed in tableland the resulting 
v(0) are shown in figure Q] With parameters constrained only by 
photometric and CO observations, the model successfully predicts 
the velocity profiles measured from the independent maser obser- 
vations. The good agreement is obtained even though the model 
profiles were not fitted to the data with any adjustable free param- 
eters. This supports the prediction given by eq.Q] as well as the 9q 
vs. Fboi relationship given by eq.|2] 
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Figure 1. Observations and modeling of the velocity profiles of three O- 
rich supergiants. The profiles are traced by various masers observations, 
shown as symbols with error bars (from Richards & Yates 1998). The model 
predictions from eq.[T]are shown with the solid lines. The terminal velocity 
■Uoo is taken from CO observations, marked by dashed lines with the la- 
uncertainty range marked by dotted lines. 



3 CORRELATIONS AMONG DYNAMICAL QUANTITIES 

Specifying the dust chemical composition determines its optical 
properties, with the relevant parameters listed in table Q] An ad- 
ditional, independent property is the dust abundance, conveniently 
parametrized with the dust cross-sectional area per gas particle 
upon grain condensation 
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Table 2. Measured and model-derived parameters for three stars with 
spatially resolved kinematic observations from Richards & Yates (1998). 
Sources are: a) HIPPARCOS database; b) Richards & Yates (1998), and 
references therein; c) from integrating observations listed in CIO (Gezari et 
al. 1993); d) L is calculated from Att D 2 F-^ \\ e) M. Rupen (private com- 
munication); f) IRAS database. 



that is, hh and rid are the number densities of hydrogen nuclei and 
dust particles, respectively, before the drift sets in (see paper I, eq. 
5). The cross-section B is a free parameter that controls the relation 
between optical depth, mass-loss rate and luminosity as (paper I, 
eq. 55) 
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where M_ 6 = M/(10" 6 M© yr -1 ), L 4 = i/(10 4 L ) and a = 
0.34Qv(T c 4 3Q7 2 *o 1 ) 1/3 , withT c3 = T c /(1000 K). While B has 
no bearing upon the shape of the radial profile v/v x , it sets the 
scale of the outflow terminal velocity, which is related to ry and 
M as (paper I, eq. 52) 
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where A = 67.5(T c 4 3 Q,*o 1 ) 1/3 X cr 2 ^ 3 . Combining eqs.[5]and 
[6] it is evident that when M is varied at a fixed L, the terminal 
velocity reaches maximum when tv = 1, i.e., at a mass loss rate 
M(w max ) oc L 3 / 4 . From this maximum, Voo decreases when M is 
either increasing or decreasing away from Af (w max ). This decrease 
reflects the role of the drift in one direction, reddening the other. 



3.1 Drift-dominated regime 

When M decreases away from M(w max ) so that tv becomes less 
than unity, the dust and gas decouple and the velocity decreases 
too. In this drift dominated regimqj, eq.|6]becomes 



T V < 1 



, ,Vt/3 , 

AM_ 6 kms 



(7) 



It is remarkable that, even though the wind is driven by radia- 
tion pressure, its velocity is independent of luminosity for small 
M (optically thin limit; see also §4.2 in paper I). This surprising 



1 For grains smaller than about 0.01 u.m, the dust drift can be negligible 
even when tv < 1; for discussion see Section CI in paper I. 
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Figure 2. Histogram of A = v oc /M 1 ^'' > (in km s — 1 ) for 36 oxygen-rich 
stars with tv < 1 outflows, with data taken from Young (1995). The mean 
(A) and its dispersion a a are listed, and the dashed line shows the Gaussian 
drawn with these parameters. 
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Figure 3. Same as figure [f] except for 60 carbon-rich stars. The data 
taken from Olofsson et al. (1993). 



result was discovered observationally by Young (1995) in a sur- 
vey of 36 nearby Mira variables with low mass-loss rates. Young 
found a clear, strong correlation between outflow velocity and 
mass-loss rate, but independent of luminosity. His correlation can 

be parametrized as Voo oc M 1 ^ 3 ' 35 , in agreement with eq. [7] if 
the observational errors for the power-law index are at least 10%, 
which is plausible. Figure [2] presents Young's data as a histogram 

• 1/3 

of Voa/M . The figure also lists the distribution mean (^4) and 
its root-mean-square scatter a a, and plots the Gaussian with these 



parameters. Figure[3]presents a similar analysis of C-rich stars with 
small optical depths, with data from Olofsson et al. (1993). In the 
case of Young's data we did not introduce any cuts since the sam- 
ple is dominated by optically thin envelopes, as can be seen from 
figure [5] On the other hand, the 63 C-rich stars in the Olofsson et 
al data include three objects with tv > 1, as is evident in figure 
[6] and these were excluded from the histogram. Each histogram 
shows a pronounced peak — in agreement with eq.|7]when the dust 
properties do not vary much within the sample. The ratio a a /{A) 
is < 25% for each sample, a fractional scatter consistent with the 
measurement errors (~10% for Vao and ~50% for M; see Ap- 
pendix for details). These strongly peaked distributions affirm the 
central role of dust drift at small mass loss rates and indicate a close 
similarity of dust properties within each sample. 



3.2 Dust-to-gas ratios 

The velocity scale A, 19 km s -1 for C-rich and 13 km s _1 for O- 
rich stars, is a fundamental property of dusty winds, derived di- 
rectly from the data. Adding assumptions regarding the grain prop- 
erties enables determination of the dust geometric cross-section per 
gas particle from 

I3A 3/2 

1/2 



13 



(8) 



where = 1.8x 10 _3 tf J /2 Q;r 1/2 T C 3 2 . The top panel of figureg] 
shows the variation with grain size of the inferred values of 6 for 
amorphous carbon and silicate grains. Remarkably, the two samples 
of different type stars produce values of 8 that agree to within 50% 
at all grain sizes. In spite of the large differences in atmospheric and 
grain properties between O- and C-rich stars, the fraction of mate- 
rial channelled into dust is such that the dust area per gas particle 
turns out to be roughly the same in both. 

The quantity closest to B that is directly determined in obser- 
vation is the dust extinction per H column density. The Galactic 
average yields Nh/tv = 1.95xl0 21 cm -2 (e.g., Sparke & Gal- 
lagher 2006), which translates to QyE> = 5. The middle panel of 
figure [4] shows the variation with grain size of the value of QyR 
inferred from the measured values of A for the C- and O-rich stars. 
Galactic interstellar dust contains a mixture of carbon and silicate 
grains of various sizes, but the dusty wind results are significantly 
below the Galactic average for both types of dust and any grain size; 
averaging the results with the MRN (Mathis, Rumpl & Nordsieck 
1977) size distribution yields the values indicated in the figure with 
dot-dashed lines. 

Another indicator of dust abundance is the dust-to-gas mass 
ratio rd g ; it is widely used even though its determination invokes 
assumptions about molecular abundances that bring additional un- 
certainty. Estimates of the Galactic average of rd g range from ~ 
0.005 (Draine 2009) to 0.01 (Barbara et al 2004). Recalling that 6 
characterizes the base of the outflow where the gas and dust have 
the same velocity (prior to the dust drift; see eq. |4), the outflow 
dust-to-gas mass ratio is 
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where p B is the density of the grain material (3.28 gem -3 for sili- 
cate dust, 2.2 gem -3 for carbon grains). The bottom panel of fig- 
ure [4] shows the variation with grain size of the value of r^g for 
dusty winds, providing yet another display of two properties noted 
above: (1) In spite of the large differences in their formation prop- 
erties, the mass fraction of carbon and silicate grains is roughly the 
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same in both types of stars, and (2) the wind dust abundance is 
significantly below the Galactic average. 

These results have significant consequences regarding the ori- 
gin of interstellar dust. Using either 8 or rd g as indicators, the dust 
abundance in winds around evolved stars is substantially lower than 
the Galactic averages irrespective of grain size or chemical compo- 
sition. The MRN average for QvB in dusty winds is more than 
a factor of 10 below the Galactic average for silicate dust, a fac- 
tor of 5 for carbon. Since the winds fail to produce the observed 
ISM values for both types of grains at all sizes, mixture averag- 
ing cannot bring the results to the Galactic values. The implica- 
tion is that dusty winds cannot be the source of all Galactic dust: 
If all dust were formed by cycling interstellar gas through stars, 
with existing grains destroyed during star formation and reformed 
around evolved stars, than the Galactic average could not exceed 
the dusty winds value. Draine (2009) discusses additional grain de- 
struction mechanisms in estimating the Galactic dust budget, and 
concludes that they could not be offset by formation in stars (see 
also Zhukovska, Gail & Trieloff 2008). Our results show that even 
without considering any processes other than cycling through stars, 
most interstellar dust is not Stardust and must have formed in the 
ISM, in strong support of Draine's conclusion. 

It should be noted that the analysis here provides a robust 
derivation of the dust abundance. The parameter A is determined 
directly from the data, and its conversion to 6 and r c \ g involves only 
minimal assumptions about the grain properties; significantly, no 
assumptions are made about any molecular abundances. The pre- 
sented results employed T c = 700 K, and the inferred dust abun- 
dance scales as r c _2 . Varying T c in its likely range, the dust abun- 
dance would increase at every a by 36% for T c = 600 K and de- 
crease by 23% if T c = 800 K; such variations have little effect on 
our conclusions. In addition, there is no need to consider global bal- 
ance of dust forming and destruction processes in the ISM, which 
are very uncertain. Our conclusion is derived from consideration 
of individual stars and therefore it is independent of, and supports, 
Draine's arguments. 



3.3 Optically thick winds 

As M increases away from A/(u max ), the wind becomes optically 
thick and reddening degrades the efficiency of the radiation pres- 
sure force, thus the velocity again decreases. The wind velocity is 
no longer independent of luminosity. Instead, from equations[6]and 
[5] it follows that when tv > 1, the final velocity is proportional to 

L 1 / 2 /M ^ 3 . Therefore, the ratio Voo/M ^ is no longer constant, 
instead it decreases as t^ 1 ^ 2 oc (M 4 3 /L) _1//2 from its value in 
the optically thin regime. Figures [5] and [6] show the comparison of 
model predictions with observations, by plotting the data and the 
relationship in eq.[6]in the u 00 /M 1/ ' 3 — M* 3 /L plane. Since the 
single free parameter A is taken from the histograms in figures [2] 
and[3] the agreement between model and data displayed in figures\5\ 
and\6\is obtained without adjusting any parameters; the coefficient 
A controls both the constant asymptotic value of Voo / M in the 
small M /L regime as well as the value of M / L at which 

• 1/3 

Vao/M begins to decrease, thus this agreement represents an- 
other test of the model. This test would be much stronger if the 
samples contained more stars with very large values of ry (i.e., 
M 4,/3 / L) so that the transitions from optically thin to thick regimes 
were better defined. We note that the samples are dominated by 
optically thin envelopes, thus our above estimates of the mean A 
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Figure 4. Dust abundance in winds around evolved stars and in the Galactic 
interstellar medium (ISM). Top panel: The dust geometric cross-section per 
gas particle inferred for dusty winds from eq.[8]as a function of grain size. 
Results for amorphous carbon grains are shown in blue solid line, for sili- 
cate in red. The inferred values for the MRN grain size distribution (Mathis, 
Rumpl & Nordsieck 1977) are shown with horizontal dot-dashed lines. Mid 
panel: Analogous to the top panel, except that it shows the product Qv$. 
whose value in the interstellar medium is directly constrained by observa- 
tions. The ISM value Qv'fi = 5 is marked with the dashed horizontal line. 
Bottom panel: Variation with grain size of the dust-to-gas mass ratio (eq.|5J 
in dusty winds. The ISM range is 0.5-1%, marked by the dashed horizontal 
lines. 

are not appreciably biased by the few stars with moderate optical 
depths. 



4 PHYSICAL DOMAIN OF DUSTY WINDS 

The relation MVoo L/c has often been invoked as the momen- 
tum conservation bound on radiatively driven mass loss rates, even 
though the mistake in this application when tv > 1 has been 
pointed out repeatedly (e.g. Ivezic & Elitzur 1995). Instead, the 
proper form of momentum conservation is Mvoo = rpL/c, where 
tf is the flux-averaged optical depth. Since rp can exceed unity for 
plausible values of tv (see paper I, eq. 57 and figure 6), momen- 
tum conservation does not impose a meaningful constraint on dusty 
winds. Instead, the constraints come from force considerations — 
the radiative outward force must exceed everywhere the gravita- 
tional pull of the star with mass M, This condition breaks down at 
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Figure 5. A comparison of the model prediction for Voa (in km s — 1 ) with 
observations of O-stars. The line plots the relationship in eq. [6] with the 
A value determined from the small-Af limit (figure|2j- The data are from 
Young (1995, dots) and Richards & Yates (1998, triangles). 



the two ends of the mass-loss-rate range, where the outward force 
is reduced for two different reasons: at very low M the gas-dust 
momentum coupling weakens thus reducing the force on the gas, 
and at very high M the coupling to the radiative force diminishes 
because of the enhanced reddening that follows increased obscura- 
tion. In paper I (see sec. 5, in particular figure 7) we derive the re- 
sulting phase-space boundaries with the aid of the appropriate scal- 
ing variables. Here we reproduce the results in terms of the system 
physical parameters. 

Figure [7] shows the lower bound on M, arising from the re- 
quirement that radiation pressure on the dust should provide suffi- 
cient force on the gas to generate a net acceleration at the base of 
the outflow. This liftoff condition sets a lower bound on M, pro- 
portional to M 2 1 L (paper I, eq. 69); below this minimal M m i n 
the grains are ejected without dragging the gas with them because 
the density is too low for efficient gas-dust coupling. As a condi- 
tion on the wind initiation, this bound is the most uncertain part of 
our solution. The result M m i u oc M 2 /L is reasonably secure (a 
similar relation was noted by Habing et al 1994), but the propor- 
tionality constant can be determined accurately only from a more 
complete formulation that handles properly grain growth and the 
wind launching mechanism. 

The bound shown in figure [8] reflects the weakening of radia- 
tive coupling in optically thick winds because of the radiation red- 
dening. For any Eddington ratio (the horizontal axis), the figure 
shows the upper limit on the wind optical depth tv (oc M A ^ 3 /L; 
see eq.|5](, determined from the full numerical solution of the self- 
similar problem; the vertical axis units were chosen to vary linearly 



_ 1 1 1 


I I J 


Model prediction 


for A=19 " 






• » • • 

• * 

* 

• 

- 


• 

• • 

• 

* •*\ 

•\ 

\* 


63 C stars from 




fi nfccnn ^ f nl i 


1 QQ ^ V 






" 0.001 0.01 0.1 

I I I 


1 10" 

I I 


0.001 0.01 0.1 


1 10 1 c 


My 6 yi_ 4 





00 



Figure 6. Same as figure[5] except for C-stars. The A value is from figure 
[5] the data are from Olofsson et al. (1993). 
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Figure 7. Winds driven by radiation pressure on dust should fall above the 
liftoff bounds shown by lines. The observations for stars with mass esti- 
mates available from SIMBAD are shown by symbols. These estimates may 
have large uncertainties due to heterogeneous sources. 



with M. Although some stars have AIvoo > L/c, none violate the 
upper bound set by proper solution of the dusty wind problem. 



5 THE WIND IR EMISSION 

Given the grain properties, the dusty wind problem requires three 
independent input parameters, which can be selected as the initial 
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Figure 8. Winds driven by radiation pressure on dust should fall below the 
reddening bounds shown by lines. The observations are shown by symbols. 



velocity, the Eddington ratio and the overall optical depth (paper 
I). However, the dependence on the first two is limited to the im- 
mediate vicinity of the boundary of the phase-space for physical 
solutions; away from that boundary, dusty winds are described by 
a set of similarity functions of the single independent variable tv. 
Figures [7J and [8] show that most stars are located well inside the 
allowed region of phase space and thus should be well described 
by tv alone. Therefore, we should be able to characterize every 
property of the wind IR emission (SED, colors, etc.) with tv only. 

We employed the code DUSTY (Ivezic, Nenkova & Elitzur 
1999) to compute the SED of each of the three oxygen-rich stars 
whose velocity profiles, as well as the approximate analytic solu- 
tion of eq.[TJ are shown in figureQ] The DUSTY calculations were 
done using the code's option "analytic radiatively driven wind", 
which computes the density profile from the same approximate ana- 
lytic solution (for details, see the DUSTY manual). Figure [9] shows 
the model SEDs, fitted to the data with the single free parameter 
Ty. The same models that successfully explain the velocity profiles 
measured from maser observations also produce satisfactory fits to 
SED of each star. The best-fit values for tv are listed in table|2] An- 
other example where SED fits provided successful predictions for 
the spatially resolved velocity profile is W Hya. Zubko & Elitzur 
(2000) fitted simultaneously the SED and the velocity profile de- 
duced from observations of the CO thermal emission and various 
masers. Significantly, neither L nor M or dust-to-gas ratio were in- 
put parameters in that model. Instead, these quantities were derived 
from general self-similarity relations after the SED fitting results 
were supplemented by the distance and velocity scales. 

Successful SED fitting is not limited to stars with silicate dust. 
We have shown elsewhere that the SED of the very dusty carbon- 
rich star IRC+10216 can be fitted as a function of pulsation phase 
by simply varying optical depth (Ivezic & Elitzur 1996). 

5.1 A Comment on Estimating Mass-loss Rate from IR 
Observations 

The infrared emission from dust is related to the gas mass-loss rate 
and may offer observationally convenient way to estimate the latter. 
The literature is abundant with various proposed expressions that 
relate IR observables and mass-loss rate (e.g. van Loon 2008, and 
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Figure 9. The observed SEDs (symbols) and model fits (lines) for the three 
O-rich supergiants whose velocity profiles are shown in figure [T] Data are 
taken from Gezari et al. (1993). The model parameters are listed in Table 2. 



references therein; see also van der Veen & Rugers 1989). How- 
ever, in addition to heterogeneous data sets and methods (see Ap- 
pendix for a summary), it is often unclear what assumptions are 
made, and what variables are considered to be independent. A key 
point of our analysis is that there are only two independent rela- 
tions between various relevant quantities, and they reflect energy 
and momentum conservation (eqs.[5]and[6j. Given these two rela- 
tionships, only two quantities can be derived, and all others have 
to be assumed, or measured. Thanks to its scaling properties, the 
problem of connecting gas mass-loss rate to IR observables can be 
decoupled into two independent steps: determination of tv from IR 
observables, and relating tv to dynamical quantities, including gas 
mass-loss rate. 

The most robust and accurate method for estimating tv from 
observations is fitting of a well-sampled SED. When the data are 
sparse, various IR colors can also be used, albeit with deterio- 
rating accuracy. Examples of such relationships are show in fig- 
ureQTJ] Given an estimate of tv, the model uniquely predicts corre- 
lations between tv and various combinations of dynamical quan- 
tities formed using v^, M, L, and fi. Examples of such rela- 
tions are shown in figure QT] (based on eqs. [5] and [6] as well as 
Mvoo = tfL/c, with the latter derived from, and not independent 
of, the first two expressions). The most appropriate expression to 
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Figure 10. The dependence of various IR colors on optical depth ry, for 
two types of dust. Predictions for other bandpasses, or for different dust 
properties, can be easily generated with the DUSTY code. 
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Figure 11. Model predictions for correlations between various dynamical 
quantities and optical depth tv (only two relations are independent). Top 
panel: Plot of the relation in eq. [6] Mid panel: Ratio of mechanical to ra- 
diative momentum. Bottom panel: Plot of the relation in eq.|5] In all panels 
Voa is in km s — , M in Mq yr _1 . The gas-to-dust mass ratio is defined as 
r g d = i/rdg (see eq.|5J. 



use for estimating gas mass-loss rate depends on the available data. 
For example, in case of the LMC studies, where outflow velocity 
is available for a much smaller number of stars than photometry, 
eq. [5] provides a superior approach because it involves only SED 
fitting (distance to the LMC can be considered well constrained in 
this context): for fixed values of a and 022, only tv and L4 are 
required to estimate M. The values for a and 022 ( see definitions 
after eq.[5j can be taken as equal to the Galactic case (Table 1), or 
could be determined as in Section [3T| using a sample of LMC stars 
with outflow velocity and mass-loss rate measurements. 



6 DISCUSSION 

We have presented here a summary of the Elitzur & Ivezic (2001) 
similarity solution in terms more suitable for direct comparison 
with observations. The correlation of and M in optically thin 
winds (eq. |7]l, first discovered observationally by Young (1995), 
emerges as a fundamental property of radiatively driven dusty 
winds and a powerful tool in the analysis of their data. Our own 
analysis in II uncovers two major new results: (1) The peaked 
distributions in figures [2] and [3] indicate there is little variation in 
dust properties among O-rich and C-rich stars, and (2) the dust 
cross section per gas particle (6) and the dust-to-gas mass ration 
(rd g ) are essentially the same for the two classes. It is hard to think 
of another method that could produce these conclusions with a sim- 



ilar level of confidence. Both of these results present a major chal- 
lenge for theoretical studies of dust formation. 

The relatively small scatter around Young's correlation sug- 
gests also a new method for determining mass-loss rates in opti- 
cally thin outflows from the relation 

M= in" 6 (^i) 3 M yr-\ (10) 

where Voo is in km s _1 and A is from £|3 . 1 1 Here M is deter- 
mined from the single measurement of Voo, without any assump- 
tions about dust abundance. In addition, there is no distance depen- 
dence, nor any need for complex modeling. The scatter in this rela- 
tion is expected to be ~ 50-70%, no worse than any other method 
for determining mass loss rates (see Appendix A). The only restric- 
tion is that the wind optical depth at visual must be less than unity, 
a condition that is easy to verify observationally. Of course, the 
absolute uncertainty of M scale is inherited from Young's (1995) 
calibration of his CO measurements, via the values of the A param- 
eter. 

In the broader context of Galactic ISM dust, our analysis rep- 
resents independent support for Draine's (2009) conclusion that 
most ISM dust particles were formed in situ, rather than produced 
in AGB winds. If verified, this conclusion would have important 
consequences for our understanding of galaxy evolution. 

Our model captures all the scaling relationships among the 
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main observable quantities that follow from energy and momentum 
conservation. These correlations can be used to estimate quantities 
whose measurements might not be available, notably the mass-loss 
rate, and for modeling the impact of AGB stars on their host galax- 
ies. An example is the recent extensive modeling study by Marigo 
et al. (2008). Our results are also suitable for analysis of massive 
data sets such as the recent SAGE survey of the LMC (Meixner 
et al. 2006). We emphasize that our derived values of B and ra g 

are directly proportional to the u 00 /M 1/ ' 3 ratio in the optically thin 
domain. If one wished to compare rd g for two populations of stars, 
say from the Galaxy and the LMC, a robust method is to compare 
the Vaa /M distributions for samples verified to be optically thin 
using spectral energy distribution. 

Despite the successful confrontation of our model with obser- 
vations, reality is more complex. The model assumes steady state 
and spherical geometry, but there is evidence that in some stars 
mass-loss rate is not steady (see, e.g., Marengo, Ivezic & Knapp 
2001, and references therein) and that some contain an additional 
bipolar component (see, e.g., Vinkovic et al 2004, and references 
therein). In addition, the physical and chemical details at the base 
of the acceleration region are not included in our model. It would be 
prudent to compare predictions for the velocity spatial profile and 
relationships among dynamical and spectral quantities produced by 
more elaborate models, e.g., Fleischer, Winters & Sedlmayr 1999, 
Hofner 1999 and Dorfi et al 2001. Such a comparison would help 
identify which features in these complex models are simply direct 
consequences of energy and momentum conservation, and which 
are unique to detailed modeling of various physical and chemical 
effects. 
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APPENDIX A: SUMMARY OF OBSERVATIONS 

Here we summarize the main measuring methods for most relevant 
observables, with emphasis on their accuracy and scaling with dis- 
tance. 

Al Dynamical Quantities 

Observational methods for determining mass-loss rate and outflow 
velocity were analyzed and compared by van der Veen & Rugers 
(1989, hereafter vdVR). More detailed discussions are presented 
by Habing (1996), Olofsson (1996, 1997) and Wallerstein & Knapp 
(1998). There are three widely employed methods: 

• The strength and shape of thermal CO line profiles contain in- 
formation about the outflow velocity and the total amount of CO in 
the circumstellar shell. The shape and width of the line profile con- 
strains the outflow velocity in an almost model-independent way. 
With the current observational capabilities and moderate signal-to- 
noise ratios, the outflow velocity can be constrained to better than 
10%. Most of the CO emission comes from the outer parts of the 
envelope, and thus the measured velocity usually corresponds to 
the final outflow velocity. 

The relationship between the implied CO mass and the directly 
observed quantities is a complex model-dependent function. The 
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expressions most often used in data analysis were derived by Knapp 
& Morris (1985). The transformation from the CO mass to gas 
mass-loss rate requires further assumptions about the outflow, its 
geometry, and the CO-to-gas ratio. Assuming a steady-state out- 
flow and spherical geometry, vdVR derive expressions that can be 
used to estimate gas mass-loss rate with an uncertainty of about 
a factor 2-5. More accurate estimates can be obtained by detailed 
modeling of individual sources, but probably not significantly bet- 
ter than a factor of 2. We note that the mass-loss rate estimate scales 
with D 2 , where D is the source distance. 

• The OH( 1612MHz) maser line profile contains information 
about the outflow velocity and the amount of OH molecules. The 
line profile usually has two well-defined peaks whose velocity sep- 
aration typically constrains the outflow velocity to better than 10%. 
This estimate is less model-dependent than the estimate based on 
the CO line profile, but unfortunately can only be used for oxygen- 
rich (O) stars. 

Baud & Habing (1983) proposed a simple model-dependent re- 
lation that can be used to estimate the mass-loss rate from the peak 
flux of OH emission. As pointed out by vdVR, the uncertainty of 
this estimate can be as large as a factor of 5 due to strong tem- 
poral variations in the OH emission strength. Intrinsic accuracy of 
this method is probably not much better than a factor of 2 due to 
a number of assumptions made, and due to uncertain values of the 
OH-to-gas ratio. We note that the mass-loss rate estimate proposed 
by Baud & Habing scales with D. 

• The third widely employed method for estimating mass-loss 
rate is based on infrared observations of dust emission. The pro- 
posed expressions (Herman et al. 1986, Jura 1987) relate the ob- 
served emission to the optical depth, which in turn is assumed to 
encode information about the gas mass-loss rate and outflow veloc- 
ity. Similarly to the above two methods that rely on the assumptions 
about the CO-to-gas and OH-to-gas ratios, the infrared-based mass- 
loss rate estimate is greatly affected by uncertain dust-to-gas ratio. 
Depending on numerous additional assumptions employed by vari- 
ous authors, the distance dependence of the infrared mass-loss rate 
estimates varies from proportional to D 2 , to no dependence at all. 

Various observables that relate infrared emission to optical depth 
have been also been widely utilized. Both Herman et al. and Jura 
assume that the 60 \im dust emission is optically thin, leading to 
t oc Feo/Fboi, where F$o is the IRAS flux at 60 |J.m and Fboi is 
the bolometric flux. Van der Veen and Rugers (vdVR) relate the 
optical depth to the F25 / F12 flux ratio, where F12 and F25 are the 
IRAS fluxes at 12 and 25 |i.m. 

One of the most widely used prescriptions for determining mass- 
loss rate from IR observations was proposed by Jura (1987). He 
assumed that the dust mass-loss rate is proportional to dust out- 
flow velocity and dust optical depth, and that dust outflow veloc- 
ity is proportional to gas outflow velocity. Dust optical depth is 
assumed to be linearly proportional to the far-IR flux (IRAS 60 
\im bandpass) because both the dust optical depth and the stellar 
contribution to the overall flux are very small at such long wave- 
lengths. Given the importance of dust drift in optically thin regime 
discussed in i)3.1| the success of Jura's formula is very surprising: 
for small optical depths the ratio of dust and gas velocities is much 
larger than unity, which should lead to significant underestimate of 
mass-loss rate. 

It turns out that there are two effects that offset each other, and 
the Jura's expression for M coincidentally produces correct val- 
ues over a large dynamic range of mass-loss rate. In addition to 
neglected dust drift, the assumption about far-IR flux being domi- 
nated by dust emission also breaks down in optically thin regime. 
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Figure Al. The ratio of mass-loss rate given by our model and that derived 
using Jura's formula is shown by the solid line, as a function of the former. 
The ratio of dust and gas velocities is shown by the dashed line, and the 
ratio of 60 urn dust emission and total flux (i.e. including the stellar contri- 
bution) is shown by the dot-dashed line. As a result of the opposite trends, 
Jura's formula is coincidentally correct to within a factor of ~2 over a large 
range of mass-loss rate. The horizontal lines are added to guide the eye and 
represent the median M ratio in the range 10 — 8 to 10 -4 Mq yr _1 (dotted 
line) and 10~ 7 to 10~ 4 Mq yr _1 (thin solid line). 

These effects are quantitatively illustrated in figure [ATI at small 
mass-loss rates the increase of the ratio of dust and gas velocities 
is well matched by the decrease of the dust emission contribution 
to the 60 tun flux. As a result, Jura's mass-loss rate estimate never 
deviates by more than a factor of 2 from the median mass-loss rate 
ratio shown by the dotted line in figure I Al I (this median ratio is 
~2.5 and represents an overall systematic offset of two mass-loss 
rate scales). Hence, the success of Jura's formula, which does not 
incorporate dust drift, is not an argument against our model which 
includes dust drift. 



A2 Photometric Quantities 

The dusty envelope absorbs the stellar radiation and reradiates it at 
longer 

wavelengths, thus making the infrared emission the most im- 
portant part of the SED for model testing. The largest catalog of 
infrared observations is compiled by Gezari et al. (1993). Individ- 
ual fluxes may often be more accurate than 10%, but due to the 
inhomogeneous nature of the catalog, the mean overall accuracy is 
probably lower. The wavelength coverage greatly varies among the 
sources and often is based only on the IRAS catalog. Recent large 
scale sensitive digital surveys (e.g. infrared 2MASS, Skrutskie et 
al. 1997; optical SDSS, York et al. 2000) are bound to significantly 
improve the availability of accurate multi-wavelength photometry. 

Given the photometric data accurate to within 10%, the bolo- 
metric flux could be determined with the same accuracy, at least 
in principle. In practice, however, the wavelength coverage can be 
sparse and this shortcoming can lead to severe errors, unless the 
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shape of SED is known a priori. Another difficulty is the variability 
of AGB stars which can also contribute significantly to bolometric 
flux errors. Given all these uncertainties, the bolometric flux can be 
determined to better than 20-30% only for a small number (< 100) 
of well observed stars. 

A3 Distance 

Several methods are employed to estimate distances to AGB stars. 
The simplest one assumes that all AGB stars have luminosity of 
10 4 Lq, and determines distance using bolometric flux. Apart from 
a bias in this estimate (the median AGB luminosity is at least a fac- 
tor of 2 smaller, see e.g. Habing et al. 1985, and Knauer, Ivezic 
& Knapp 2001), its intrinsic accuracy cannot be better than about 
20-30% due to the finite width of the AGB luminosity function 
(Jackson, Ivezic & Knapp 2002). Other methods that have been 
frequently used to estimate distance include period-luminosity re- 
lations (e.g., Whitelock, Marang & Feast 2000), assumption that 
the absolute K-band magnitude is the same for all stars, and kine- 
matic estimates based on radial systemic velocities. The distance 
errors associated with these methods are hard to characterize and 
sometimes do not even have a Gaussian distribution (e.g. kinematic 
distances); they may be more accurate than a factor of two, but 
probably not better than 20-30%. 

The distance estimates to nearby AGB stars have been recently 
greatly improved with the data obtained by the HIPPARCOS satel- 
lite. The accuracy of HIPPARCOS distances varies from a star to 
star, but, nevertheless, there are now hundreds of AGB stars with 
distance estimates better than 10%. 

A4 Bolometric Luminosity 

The bolometric luminosity can be determined using the bolomet- 
ric flux and distance estimates. Assuming HIPPARCOS distances 
and stars with good photometric coverage, the luminosity can de- 
termined to within 10-20%. In more typical cases, its uncertainty 
is closer to 50% (e.g. Knauer, Ivezic & Knapp 2001), and without 
a Hipparcos parallax it can be as large as a factor of 2 (Jackson, 
Ivezic & Knapp 2002). 
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